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Abstract 

We study the nonequlibrium transport through a quantum dot weakly cou- 
pled to Luttinger liquids (LL). A general current expression is derived by 
using nonequilibrium Green function method. Then a special case of the dot 
with only a single energy level is discussed. As a function of the dot's energy 
level, we find that the current as well as differential conductance is strongly 
renormalized by the interaction in the LL leads. In comparison with the sys- 
tem with Fermi liquid (FL) leads, the current is suppressed, consistent with 
the suppression of the electron tunneling density of states of the LL; and the 
outset of the resonant tunneling is shifted to higher bias voltages. Besides, 
the linear conductance obtained by Furusaki using master equation can be 
reproduced from our result. 
PACS numbers:73.63.-b,71.10.Pm, 73.40.Gk, 73.23.-b, 
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Tunneling phenomena in mesoscopic systems have been intensively studied for years both 
experimentally and theoretically. Since early 1990s, the tremendous progress in nanofab- 
rication and the discovery of novel one-dimensional (ID) materials, especially the carbon 
nanotubes, have led to a revived interest in studying the tunneling in systems including 
ID conductors. § [0 

It is well-known that at low temperatures the two-dimensional (2D) and three- 
dimensional (3D) interacting electron systems can be described by Landau's Fermi liquid 
(FL) theory. However, the interacting electron systems in one dimension (ID) is quite differ- 
ent, it behaves as Luttinger liquid (LL), ^ for that theory predicts a series of unique 
characters, e.g., the absence of Landau quasiparticles, spin-charge separation, suppression 
of the electron tunneling density of states (DOS), anomalous power laws for transport co- 
efficients with interaction-dependent exponents, etc. p[ Naturally, one may ask whether 
the (resonant) tunneling between LLs is different from that of FLs. 

Recently, Furusaki reported a theoretical study of the resonant tunneling through a quan- 
tum dot (QD) weakly coupled to two LL leads (hereafter as LL-QD-LL). ||^ Assuming that 
the tunneling is an incoherent sequential tunneling, the author derived the linear conduc- 
tance using the master equation approach. Similar studies have been done by other authors 
earlier. [|T^] |TT| In this paper, we present an investigation for a LL-QD-LL system. 
Differently, we consider the tunneling as a coherent process and study the nonequilibrium 
current through the system with a finite bias voltage applied on the two LL leads. By 



using nonequilibrium Green function technique, []T3| a general current formula is derived. 
We find that both the current and differential conductance are strongly renormalized by 
the interaction in the LL leads. In current-voltage characteristics, the current is first sup- 
pressed, consistent with the suppression of the electron tunneling density of states; and the 
outset of the resonant tunneling is found to be shifted to higher bias voltages, comparing to 
the system FL-QD-FL. Finally, the linear conductance obtained by Furusaki using master 
equation can be reproduced from our result. 

The Hamiltonian of the system can be split into three parts: H = Hieads + + Ht, 
where Hi^ads = Hi + denotes the left and right LL leads in its standard form (see 
0), = Ylin^ndndn IS the Hamiltoniau of the QD, with dn} the creation/annihilation 
operators of the n-th energy level in the dot (for simplicity, the intradot Coulomb interaction 
has been neglected). Ht is the tunneling Hamiltonian and can be written as 
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Ht= Yl {tn^xd+tPx + h.c), (1) 

n,\=L/R 

in which {V'a ' V'a}('^ = L/R) are the Fermi operators at the end points of the left/right lead; 
{tn,\} is the tunnehng constant. The spin indices have been suppressed since they are not 
important in the following discussion. 

Defining the mixed casual Green function Gn^x (^1,^2) = ~i {Tdn (ti) 'ipx (^'^))h Green 
functions of all the other types in the standard form, [|1^] the current from left LL lead into 
the QD under a bias voltage —V applied on the right lead can be evaluated as follows (in 
units of h = 1): 

-^^ = 4^) = -^'^.^^ KLG<At,t)] , (2) 

where Gn^iij^T') = Xlm /c'^''"!^".'" ''"1) ^m.Lfi'L (ti, t') is the Dyson equation for contour 
Green function Gn,L (t, t'), Gn,m (t, t') = —i {Tcdn (r) {'t'))h ^^^^ Green function of 

the QD, (r, r') = —i (Tcipi (t) ('''')) Hl local Green function at the end point of 

the left LL lead without the coupling to the QD, and is the contour-ordering operator. 

After using Langreth theorem of analytic continuation, and the Fourier transformation, 
the current can then be expressed as: 

^ [t:,L {G:,m (^) 9l (^) + G<^ (u) gl [uj]) t^, J . (3) 

The Dyson equation for Gn,m is Gn,ni (t^) = ^n,m9n (^) + E« 9n (t^) ^n,i (t^) G/,™ (w), 
where E„^; (c<j) = J2x=l/r^"',>'^Ix9>' irreducible self-energy and gn{oj) is the free 
Green function of the dot without the coupling between the QD and leads. Using Langreth 
theorem and noticing Y7'°- = ±i (E> — the current expression, Eq.(3), can be reduced 
to 



Neglecting the off-diagonal tunneling matrix elements and their energy-dependence, 
which is equivalent to take t*^itm,L = ^n.m |^n,L|^, then Eq.(3) reduces to 



which describes the nonequilibrium current in LL-QD-LL, and is the the central result of 
this work. Similarly, one can easily obtain the current from right LL lead to the QD, J^j, 
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just by changing the subscript in Eq.(5) L ^ R, and easily to see Jl + Jr = 0, the current 
conservation, as expected. 

As an important example, let us consider the quantum dot with only a single energy 
level, i.e., Hd = ed^d. After taking a contour integral, Eq.(5) reduces to 

In which g'f^'^ (e) and g"^'^ {e) are the lesser (greater) Green functions of the left and right 
leads without the couplings to the dot, and can be directly quoted from: p 

at'^ie) =±^-^e^^7L(e), (7) 

T V 
9r^ (^) = ^e^'^Ti? {e - eV) , 

where T is the temperature, and V the bias voltage. 1l/r{.^) is defined as 7L/i?(e) = 

I 1 2 

[j^Y^^^^^ ^ r('i/g^^^ ) — '' ^^"^^ ^ (^) Gamma function, gi/R are interaction 

parameters characterizing the left/right LL liquids, which should not be confused with the 
Green functions with the same symbols, Tij^ describes the effective level broadening of the 

I 1 2 

dot, proportional to , and A is the high-energy cutoff or a band width. P 

In the lowest order approximation, Eq.(6) reduces to 

J = eT lLie)^nie-eV) / eV\ 

7L (e) cosh (^) +^R{e- eV) cosh (^) \2T J ' ^ ^ 



Eq.(8) describes the nonequilibrium current of LL-QD-LL under a finite bias voltage. In 
the following, we present the numerical studies based on Eq.(8) for the symmetric case, i.e., 
9l = gR = g, ti = tn and Tl = Tr. 

Fig.l shows the current vs bias voltage at a fixed value ^ = 5 for different interaction 
parameters of the LL leads, g = 0.9,0.7, |, |, ^. The system with PL leads is also shown 
for comparison (see the curve with g = 1 in Fig.l). We find the following features: (1) For 
FL-QD-FL , a steep increase occurs around |^ = 5, which is consistent with the suppression 
of the electron tunneling of states of the LL. The suppression becomes stronger with the 
decrease of g (or the increase of the interactions in the leads). (2) The outset of the resonant 
tunneling varies with the change of g, from eV = e for g = 1 (i.e., the FL leads) to eV = 2e 
for ~ (i.e., the strongest interacting LL leads). A current plateau appears again for 
g ^ around eV = 2e. (3) Somehow, the curves with g ^ ^ intersect approximately at a 



point (|^ = 10, -p— = |), which may indicate the pecuhar behavior for strong interacting 
case. 

The second feature above mentioned can also be seen from Fig.2, in which the differential 
conductance vs ^ for different interaction parameters of LL, (7 = 1, i, ^, is presented. The 
location of the peak of conductance is shifted from e = eV to e = ^ with the increase of g. 

The nonequilibrium current as a function of ^ at a fixed bias is shown in Fig. 3. Compar- 
ing to the current of FL-QD-FL, these curves are strongly renormalized by the interactions 
of the LL leads: instead a mesa for FL-QD-FL, a peak occurs at e = and the half- 
width of the peak becomes narrower with the increase of the interactions of electrons in 
the LL leads. The dependence of the peak current on T and e can be obtained by esti- 
mating the current J at e = Because the asymptotic limit of the Gamma function is 
F (x) ~ x^~^e~^\/2Ti ( for |x| — > 00, and |argx| < n ~ S, S > 0), [l^ one finds that at low 
temperatures (T <^ |) the current behaves as J ~ {^Y^^ ^, whereas at high temperatures 
(T ^ 1) , J ~ ^ {jxY^^ ^- Similar scaling forms have been derived by Chamon and Wen 
by considering the incoherent sequential tunneling. [|1^ 

It should be pointed out that our calculation is valid only when the higher order contri- 
butions (> O [\tLtRf)) can be neglected. The pole of the integrand in Eq.(4) satisfies the 
equation {g^ (c<j))^^— S"" (uj) = 0. Neglecting the shift of the energy level of the QD, the lowest 
order solution is u = e + ImY7' {e). Thus, the higher order correction can be neglected only 



if 



ImE''(£+ImS''(E))-ImS''(£) 



7(£+ImS'-(£))-7(£) 
-lie) 



<^ 1. Using the asymp- 



^ 1, or equivalently, 

totic limit of the Gamma function, one finds that at low temperatures {t ^ ^^^r^j ■> 

condition reduces to ^ (^)^^^ ^ In (^) ^ 1 and ^ { ^^\k^^ ^ 1^ (^^^^) whereas 
at high temperatures {t ^ |, ^^^y^j , ^ (x)^^^ ^ ^ 1- Therefore, our evaluations validate 



only when T 7^ and e, eV, T <^ A. [T^ 

Finally, as a quick check, we can obtain the linear conductance a from Eq.(7) by simply 
taking V 0^, which gives: 



a = I (9) 

2cosh[^] 7l(£) + 7r(^)' 

This is exactly the linear conductance obtained by Furusaki, and was confirmed by 
experiment recently. 

In conclusion, we derive a general formula for nonequilibrium dc current through a quan- 
tum dot weakly coupled to two Luttinger liquids. We consider the tunneling as a coherent 
process and use nonequilibrium Green function technique, different from Chamon and Wen 
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and Furusaki. In symmetric case and for the dot with a single energy level, we find that both 
the current and the differential conductance are strongly renormalized by electron-electron 
interactions in the LL leads. We hope that further experiments may provide evidence for 
these theoretical predictions. 
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FIGURE CAPTIONS 



Fig. 1 Current, scaled by its maximum at ^ = 5, vs bias voltage for g — 0.9, 0.75, \ ,\ ,jq- £ 
is the energy level of the dot, g is the interaction parameters of the leads. Curves for 
g — intersect approximately at a point — 10, = . The case with FL 

leads {g — 1) is also shown for comparison.. 

Fig. 2 Differential conductance scaled by its maximum at |^ = 10 as a function of ^ for 
g — 1, |, The location of the maximum is shifted away from e — eV iov g — 1 to 
£ = ^ for ~ 0. 

Fig. 3 Current vs ^ for = 1, |, ^^j. The current is scaled by its maximum at bias voltage 
1^ = 10, For g < l(Luttinger liquid), a peak occurs at e — and the half width 
of the peak becomes narrower with the decrease of g. The curve with g — 1 (Fermi 
liquid) is also shown for comparison. 
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